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Abstract 

In this contribution we explore one of the many possibilities of deter- 
mining the strong coupling constant as at hadron colliders. The method 
considered is quite unique compared to other methods in that the value of 
as is determined by the "evolution rate" of the parton density functions 
rather than by the "event rate" . 



1 Talk given at the "5th International Workshop on Deep Inelastic Scattering and QCD" 
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Introduction and Motivation 



Hadron colliders will supply an increasing amount of data with the upcoming 
high luminosity TEVATRON run and the LHC project. Methods for extracting 
as and parton density functions from these data sets can therefore expect a 
steady improvement in the precision over the coming decades. 

Often it is claimed that hadron colliders cannot do such precision measure- 
ments. A few comments can be made in answer to this. First of all, a hadron 
collider measures the value of as at many different values of the (partonic) cen- 
ter of mass energy. This is in contrast to e + e~-colliders where the center of 
mass energy is fixed. In fact, at the TEVATRON the partonic center of mass 
energy useful for the ag-extraction can go as high as 1 TeV (and at the LHC this 
will increase by an order of magnitude) . Secondly, a hadron collider can make 
accurate measurements by selecting appropriate observables. For this purpose 
we select in this talk the normalized one jet inclusive transverse energy distri- 
bution. The value of as will be determined from the shape of the distribution. 
The major factor determining the shape is the fraction of quarks in the colliding 
hadrons. Due to the evolution of the parton density functions the quark frac- 
tion at moderate parton fractions decreases as the jet energy increases. This 
depletion is controlled by the strength of the strong coupling constant. 

Because we look at a normalized quantity both the theoretical (renormal- 
ization scale dependence) and experimental uncertainties (e.g. luminosity un- 
certainty) are much smaller than one would expect. 

Measurement and Methodology 

As mentioned in the introduction, the observable used is the normalized one jet 
inclusive transverse energy distribution. As an example we use the published 
run la results from the CDF collaboration fl. For the theoretical prediction 
we use the JETRAD monte carlo 0j with the cuts and jet algorithm as close as 
possible to the experimental setup. The MRSA' parton density functions 
which allows varying as, were used. The renormalization/factorization scale, 
/i, was chosen to be a constant, A, times the maximum jet transverse energy, 
Et, in the event. Both data and theory are divided by the "reference" theory 
prediction which is given by: as(Mz) = 0.120, [i = Et- To normalize the 
distribution we choose the ratio to be equal to unity at Et — 200 GeV . 

We show the as-dependence in fig. la and the scale dependence in fig lb. 
As can be seen the dependence on as is quite substantial compared to both 
the experimental and theoretical uncertainties. Also note that the leading order 
(LO) and next-to-leading order (NLO) results are quite close. The only differ- 
ence between the two predictions is the (expected) reduced scale dependence 
at NLO. The method to extract as is now quite simple: we minimize the x 2 
to fit the theory to the data in fig. 1 between 30 GeV < E T < 200 GeV by 
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Figure 1: (a) The sensitivity of the distribution to the value of as(Mz). For 
comparison the CDF data is also shown, (b) The scale dependence of both 
the leading (LO) and next-to-leading (NLO) predictions. 



varying both Q5 and A. The results and, more importantly the interpretation 
are discussed in the next section. Note that we do not consider the systematic 
uncertainties at this point. They can be easily included in the x 2_n t by building 
up the correlation matrix out of the systematic uncertainties Q. Their inclu- 
sion is better left to the experimenters. Here we want to concentrate on the 
methodology. 



Interpretation and Results 

In fig. 2a we show the results of the minimalization procedure to fit to the 
data. Both the minimum and the 1-a uncertainty ellipse is shown. The figure 
contains all the information we can extract from the data. While the central 
value is quite trivial to determine, the interpretation of the uncertainty is not. 
The perfect answer (that is no renormalization scale uncertainty) would be a 
vertical strip. The uncertainty would then simply be the width of the strip 
independent of the choice of A. However, in fixed order perturbative QCD we 
have a residual scale sensitivity due to the truncation of the series. This is 
reflected in the slope of the ellipse-axis. In fact one could argue that the slope is 
the correct measure of the theoretical uncertainty. When comparing results from 
different experiments this slope could be used to weight different experiments 
on their theoretical uncertainty. This all implies that the parameter A cannot be 
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Figure 2: (a) The result of the x 2 -fit to the data for both the LO and NLO 
predictions as the minimum and the 1-er uncertainty ellipse, (b) A comparison 
of the extracted central value of ag LO {Mz) with its one- two- and three-loop 
evolution compared to as(£r) extracted from the data using the method of ref. 



considered a fitting parameter as as is, nor is it in a direct manner related to the 
theoretical uncertainty. Note that when the data accuracy increases (e.g. the 
CDF/DO run lb data) the fit of the theory to the data will become more strained 
and the variation of A will become more constrained. This does not indicate 
that the theoretical uncertainty is decreasing. On the contrary, this means that 
the NLO prediction is becoming more and more inadequate to describe the data 
and even higher order calculations are needed. 

For the moment we use a naive procedure to quote the theoretical uncer- 
tainty. The experimental uncertainty is taken to be the width of the ellipse at 
the minimum, while the theoretical uncertainty is taken as the variation within 
l-<7 for scales between 1/3 < A < 3. The results are 

ol l s °{Mz) = 0.110 ± O.OOl(stat) ± 0.004(theory) 
a^ LO (M z ) = 0.114 ± O.OOl(stat) ±0.004(theory) . 

Alternatively, one could argue that the difference between the LO and NLO 
value of ag should be larger than the difference between the NLO and NNLO 
value of as, giving an alternative, but equal, estimate on the theoretical uncer- 
tainty of 0.004. 
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Conclusions and Outlook 



We have used the normalized one jet inclusive transverse energy distribution to 
extract as(Mz)- The fact that we used the normalized distribution reduces the 
experimental and theoretical uncertainty significantly. To improve the results 
wc need to fit the partem density functions (specifically the gluon) together with 
the value of as- This should remedy the obvious discrepancies between the data 
and theory for Et > 200 GeV (see ref. || ) . Such a measurement would not only 
give us as, but simultaneously a true NLO determination of the gluon parton 
density function. In fig. 2b wc finally show, as a cross check, the comparison 
between the as extracted in this talk and the as-values determined in each 
Ex-bin, obtained using the methods of ref. || (incorporating parton density 
functions with varying as)- As can be seen the agreement between the two 
methods is quite good. 

The new run lb data from the CDF and DO collaboration will severely test 
the NLO description of the data and higher order calculations might be needed 
to describe these results and extract the gluon parton density function and as- 
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